Critical behavior of the long-range Ising chain from the largest-cluster probability 

distribution 
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Monte Carlo simulations of the ID Ising model with ferromagnetic interactions decaying with dis- 
tance r as l/r 1+CT are performed by applying the Swendsen-Wang cluster algorithm with cumulative 
probabilities. The critical behavior in the non-classical critical regime corresponding to 0.5 < a < 1 
is derived from the finite-size scaling analysis of the largest cluster. 
PACS numbers: 05.50.+q, 64.60.Cn 



43 

o 



CO 



I 

C 

o 
o 



> 
o 



o 
o 



-a 
c 

o 
o 



X 



The calculation of the critical behavior of the Ising 
model with long-range (LR) interactions decaying with 
distance with a power law \jr d+cr is not an easy task, 
even in one dimension, where phase transition at finite 
temperature occurs for < a < 1 0^]. Exact analytical 
expressions for the critical exponents may be derived for 
a < 0.5 ||, corresponding to the classical regime, while 
for a > 0.5 only the approximate results exist. Various 
analytical and numerical approaches have been applied, 
from direct numerical calculations on finite chains [ 5} , to 
several approaches based on the renormalization group 
(RG) and scaling ]5|-[To[|. The nonlocal character of in- 
teractions reduces the efficiency of most of the standard 
approaches so that the values of critical exponents ob- 
tained by all these methods differ considerably. 

In numerous cases of phase transitions in systems with 
short-range interactions, a useful complementary tool for 
obtaining both qualitative and quantitative results is pro- 
vided by Monte Carlo (MC) simulations in combination 
with finite-size scaling. 

Such systematic studies were lacking for the LR models 
until recently. Namely, when applied to models with LR 
interactions, the standard MC approaches based either 
on Metropolis or on various cluster algorithms are par- 
ticularly time consuming, since the number of operations 
per spin-flip is proportional there to the size of the sys- 
tem. Recently, this problem was successfully resolved by 
Luijten and Blote |Tl]] who used the cumulative probabil- 
ities within the Wolff cluster algorithm In^jWhich they 
applied to the Ising and similar models p3[-|f 5| reducing 
the computing time by several orders of magnitude. 

Their very exhaustive studies concentrate on questions 
related to the mean-field (MF) regime, while very little 
or no interest has been dedicated yet to the regime of the 
non-classical critical behavior corresponding to a > 0.5. 

The purpose of this work is to extend the MC stud- 
ies of the critical behavior of the LR Ising model to the 
non-classical regime. At the same time this is a suitable 
example to examine the efficiency of using only the clus- 
ter statistics in deriving the critical properties of the LR 
model. 

The ID Ising model with LR interactions written in 
form of a special case (q = 2) of the Potts model is de- 



scribed by the Hamiltonian 
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where J > 0, Sj is a two-state Potts variable at the site 
i, 5 is the Kronecker symbol and the summation is over 
all pairs of the system. By the substitution S Si , Sj = (Si ■ 
Sj + l)/2, and Jj = J/2, where Si = ±1 and Jj denote 
the Ising spins and the interaction constant, respectively, 
one recovers the standard definition of the Ising model. 
In the mean-field regime, a < 0.5, the critical exponents 
have classical values (y = 1/cr, rj = 2 — a). We focus here 
on the region 0.5 < er < 1, where the critical exponents 
are non-trivial and known only approximately. 

When the shortcut using the cumulative probabilities 
is implemented in the MC simulations, the full extent 
of the numerical advantages is achieved when only the 
distribution of spins has to be calculated, while for the 
energy sampling the reduction in the CPU-time is less 
important. Recently the calculation by energy sampling 
was improved [ fl6| by calculating the energy in the mo- 
mentum space and applying the fast Fourier transform. 
Our intent here is to avoid the energy sampling by mak- 
ing a more complete study of the cluster statistics. To 
this purpose we use the Swendsen-Wang cluster algo- 
rithm |17|| . The implementation of cumulative probabili- 
ties to the Swendsen-Wang cluster algorithm is straight- 
forward. Each step of this iterative procedure consists 
in identifying all the clusters in a given spin configu- 
ration of the system following the rule that two par- 
ticles belong to the same cluster with the probability 
Pij = (1 - exp[—J(i,j)/T])6 SitSj , and then flipping all 
the clusters randomly. The cumulative probabilities are 
applied at the point of identifying the individual clus- 
ters along the same lines as was done |n|] for the Wolf 
single-cluster algorithm, and this reduces the number of 
required operations per a single spin-flip by a factor of 
system size. Like the Wolf algorithm, the Swendsen- 
Wang algorithm suppresses the slowing-down at critical- 
ity, and though it might be somewhat more costly re- 
garding CPU time, it gives a more complete insight into 
the cluster statistics and the related probability distribu- 
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tions, such as the distribution of cluster sizes or that of 
the largest cluster. 

The probability distributions related to clusters are 
mostly used in the description of critical properties of ge- 
ometrical transitions, like the percolation . As it fol- 
lows from the graph expansion by Kasteleyn and Fortuin 
, the very basis of the Swendsen-Wang algorithm, this 
cluster statistics may be related to the thermodynamic 
quantities of the thermal transitions as well. The only 
important difference is that in percolation one deals with 
the simple geometrical clusters, while in thermal transi- 
tions the clusters formed by the active bonds (in terms 
of the cluster algorithm) are considered. In systems with 
LR interactions, such as the present model, this differ- 
ence is obvious. 

The cluster-related distributions were much less ex- 
ploited in the case of thermal transitions, although some 
detailed studies exist (see e.g. refs [^0 21 1) for models 
with short-range interactions. More recently, there has 
been renewed interest in this subject within a somewhat 
different context [ p2p3[ ]. 

By performing lO^MC iterations, which is sufficient 
for most of data presented here, the chains of sizes up 
to a few tens of thousands of sites could be simulated 
without requiring too much numerical effort on a modest 
workstation. 

We limit our analysis here to the size distribution of 
the largest cluster, the quantity that was left out of reach 
within the earlier single-cluster approach. 

The probability that the largest cluster will be of the 
size I is defined as 
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where N max {l) is the total number of occurrences of the 
largest cluster of size I during umcs MC swaps of the 
system. 

Several thermodynamic quantities of interest may be 
expressed through the corresponding moments defined as 



< r >- 



The first moment 
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gives the mean size of the largest cluster. Below the crit- 
ical temperature it describes the order parameter, which 
is usually defined as the average of the largest compo- 
nent, i.e. M = (q < max{m a } > — l)/(q — 1), where m a 
denotes the fraction of spins of the system in the state a. 

A quantity of interest to consider is the Binder's 
fourth-order cumulant [E4j ratio for the largest cluster, 



defined in the same way as for the standard order pa- 
rameter. Differences between the two corresponding dis- 
tribution functions, however, result in different shapes of 
the same ratio defined for the two quantities. As may 
be seen in Fig. 1, the ratio Rl(T) at high temperatures 
and in the thermodynamic limit tends to 1 and not to 
3, while the corresponding probability distribution is ex- 
pected to be non-Gaussian, as found already in the case 
of percolation |2j| and the short-range Ising model . 
However, as well as in the mentioned cases, these dif- 
ferences should not change the basic scaling properties 
in the critical regime. At low temperatures and around 
T c the behavior of the two ratios coincides. Like in the 
case of the standard order parameter, there is a com- 
mon crossing point of the curves Rl{T) of different sizes 
where the ratio (j^) almost does not depend on L and 
which can be identified as the critical temperature. In 
Fig. 1, Rl{T) is shown in the case a = 0.9 and for sizes 
varying from 1000 to 20000. By taking into account sizes 
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FIG. 1. Ratio Rl{T) for a = 0.9 and sizes L varying from 
1000 to 20000. 



within the range 1000 < L < 20000, the determination of 
this crossing point can be done with precision of 4 digits, 
which gives the estimated error of 0.2%. 

TABLE I. Comparison of results for the critical tempera- 
ture T c to earlier results by Nagle and Bonner 0, Monroe et 
al. §, FRS §, Cannas and Magalhaes @ and Monroe @. 



this work 



H 



2t 



0.6 1.770 ±.003 1.766 1.7885 1.7718 

0.7 1.463 ±.003 1.458 1.491 1.4635 1.28 1.466850 

0.8 1.2155 ±.002 1.212 1.2585 1.2150 - 

0.9 1.001 ±.001 1.0015 1.058 1.0027 0.77 1.018845 



Rl(T) = 



< I 2 > 2 



(5) 



2 



In table I are presented the critical temperatures for 
several values of ct in the non-classical regime 0.5 < a < 1 
incremented by 0.1, obtained as common crossings of 
Rl(T). For comparison, we quote the earlier results ob- 
tained by numerical calculations [Q, coherent anomaly 
approach Pjgq] , finite-range scaling (FRS) g, and real- 
space RG |lOfl. When quoting ref. || we mention only 
one of the two given sets of results for T c , which fits bet- 
ter to the later work |^6| performed with a significantly 
improved precision, but, unfortunately, available only for 
two values of the range considered here. 

Near T c (i.e. in the regime L < £), the moments (||) 
and their ratio (Q) should obey finite-size scaling. Thus 



Rl(T) = L x f{L x ' v t) 



(G) 



where r = (T — T c )/T c is the reduced temperature. Since 
Roc(T c ) is finite, x = 0. 

It turns out that the critical exponent v may be esti- 
mated with a rather good precision by a simple tuning 
of the unknown exponent \jv in the Eq. (^|) until the 
curves collapse in the vicinity of the crossing point. 

In Fig. 2 we present the curves from Fig. 1, collapsed 
after taking l/v = 0.40. 
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FIG. 2. Ratio Rl(T) for a = 0.9 versus reduced temper- 
ature r scaled with L 1 ^ 



The results for v -1 obtained by collapsing fits for other 
considered values in the interval 0.5 < a < 1 are given in 
table y. In order to reduce the degree of arbitrariness of 
such fits we require the curves to be indistinguishable in 
the region of temperatures around T c , up to the tempera- 
ture of the maximum of R^{T), T max . This temperature 
may be related to the finite-size crossover temperature 
up to which the finite-size scaling equation (|6j) should 
hold. In Fig. one can observe that the maximum of 
Rl(T) shrinks and shifts towards T c with increasing L. 
The rough fit of the difference T max —T c to the power-law 



TABLE II. Results for u 1 compared with earlier results 
of Nagle and Bonner Q, Fisher et al. ||, Monroe et al. Q, 
FRS U and Cannas and Magalaes [jToj. For the meaning of 
labels a and b, see the text. 



a this work 



ja flb 



fla flb p| P 



0.6 0.50 ± 0.02 0.504 0.540 0.5098 0.441 0.471 0.501 

0.7 0.51 ±0.02 0.542 0.552 0.4856 0.478 0.461 0.518 0.376 

0.8 0.47 ± 0.02 0.533 0.538 0.4461 0.505 0.446 0.483 - 

0.9 0.40 ± 0.01 0.529 0.526 0.4032 0.542 0.421 0.405 0.256 



form (1/L)^ gives the values of tfi which roughly corre- 
spond to l/v for a given ct, as expected for the finite-size 
crossover exponent p^ ]. The interval of tuning in which 
the curves remain indistinguishable determines the er- 
ror margins of the obtained values for the exponent v , 
which increase with decreasing ct but do not exceed 4% 
in the considered range of cr's. 

Alternative calculations of the exponent, e.g. by using 
the derivatives of moments, turn out to be less advan- 
tageous. (The purely numerical derivative would require 
a much higher numerical precision and consequently far 
more extensive runs; the derivative obtained by his- 
togram interpolations or by mixed moments would re- 
quire an energy sampling and an increased numerical 
precision in addition.) We also mention at this point 
that the calculations of the Binder's cumulant ratio may 
equally be performed for the standard order parameter 
using the Wolf cluster algorithm. One advantage of the 
approach which uses the Swendsen-Wang cluster algo- 
rithm is that it has similar precision and efficiency both 
in low- and high-temperature regimes, which improves 
the scaling fit. 

The degree of precision of the obtained data permits 
the comparison to the results by other approximate meth- 
ods. These results, which are quite few, are quoted in ta- 
ble H in chronological order and correspond to numerical 
calculations on finite chains using Pade series extrapo- 
lations [Q, e-expansion near 2a — d ||, the coherent 
anomaly method ||], finite-range scaling ||, and real- 
space RG . In refs j3|] and Q only the numerical data 
for the exponents 7 and (3 are available. The exponent v 
may be expressed in two ways, which, due to the approxi- 
mate character of the results, do not give identical values. 
By "a" are labelled the values obtained from the scaling 
relation v = 7/(2 — 77) and the exact analytical expression 
r\ = 2 — a expected to hold also in the whole non-classical 
regime. By "b" are labelled the values obtained from the 
scaling relation involving both exponents, v = 7 + 2(3. 
The quoted results of ref. [|| are the ones obtained by 
the expansion in Act = a — 0.5 up to the second order. 
The expansion in (1 — a) || is not quoted in the table. 
It holds very close to ct = 1 and we may quote only the 
value for ct = 0.9 which is equal to v' 1 = 0.447. The re- 
sults of ref. |5| are also expected to be more reliable closer 
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to a = 0.5. Other approaches cover the whole range of 
a evenly, but their results are rather different and their 
accuracy as well. Our MC results agree the best with the 
FRS results. 

In ref. jl(|, which deals with the Potts model with 
arbitrary number of states q, we found intriguing the 
conjecture that the exponent v would remain the same 
for other values of q (when the transition is of the second- 
order), which is in contrast to earlier results of FRS ||. A 
MC approach can thus be used to provide an independent 
calculation. (It also has an advantage over both used 
approaches that it is able to distinguish the first- from the 
second-order transition j2^|jl|). Our first calculations 
p9| , performed along the same lines as above, give for 
the case q = 3, a = 0.9 the value l/v = 0.48 ± 0.01, 
which would rather be in favor of the FRS results. 

In conclusion, we have applied the MC simulations to 
the two-state Potts model with LR interactions on chains 
up to 20000 sites by using the Swendsen-Wang algorithm 
to calculate the cluster statistics, in particular the size 
distribution of the largest cluster. The approach was 
used as an alternative calculation of the critical behavior 
in this model in the non-classical regime 0.5 < a < 1. 
We have shown that the scaling analysis of the Binders 
4th-order cumulant ratio of the largest-cluster size gives 
the critical temperature and the critical exponent v with 
reasonable accuracy (better then 0.2% and 4%, respec- 
tively). Although the overlapping fit is not expected to 
be very precise, nor were the simulations pushed to their 
extreme, the estimated errors are significantly smaller 
than the differences between the results obtained by dif- 
ferent approaches. The obtained results, presented in 
tables I and II, are closest to the values obtained ear- 
lier by the FRS. While the present work deals only with 
the Binder's 4th-order cumulant ratio, a wider analysis, 
which includes other quantities that may be derived from 
the cluster statistics, should be performed in future. We 
believe that the approach through cluster statistics might 
be useful to study other cases of the Potts model and 
other models with discrete symmetries, where in the case 
of LR interactions, the alternative methods are very re- 
stricted. 
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